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ABSTRACT 


Applications  have  been  made  of  the  localized 
perturbation  technique  for  conjugated  systems  within  the 
framework  of  Htickel  theory  by  using  the  inverse  of 
cyclic  matrices  in  closed  analytical  form  in  terms  of 
Chebyshev  polynomials  by  the  method  of  Lowdin,  Pauncz 
and  de  Heer.  It  is  shown  that  by  this  approach,  one 
can  get  the  exact  energy  levels  and  eigenfunctions  of 
perturbed  conjugated  systems  by  a  simple  non-iterative 
method  for  simple  perturbations  and  by  a  second  order 
iteration  process  for  complicated  perturbations  which 
reduces  the  numerical  "size"  of  the  problem  considerably. 
Possible  applications  of  this  technique  in  other  related 
fields  are  discussed. 
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INTRODUCTION 

In  this  technical  note,  we  shall  apply  some  of  the  methods  developed 
by  Ldwdin  in  sections  VI  and  VII  of  the  "Studies  in  Perturbation  Theory" 
in  an  investigation  of  some  properties  of  perturbed  simple  conjugated 
systems.  These  applications  will  illustrate  the  power  and  elegance  of  this 
new  approach.  Most  of  the  results  contained  herein  are  not  entirely  new  but 
the  methods  used  have  been  recently  developed  amd  show  promise  of  being 
particularly  powerful  when  applied  to  more  complex  systems  than  those 
described  here.  Thus  the  emphasis  shall  be  more  on  the  techniques  than  on 
the  results. 


P.O.  Lbwdin,  Technical  Notes  Nos.  64  and  65  from  the  Uppsala 
Quantum  Chemistry  Group. 


We  shall  be  working  within  the  framework  of  the  simple  Hflckel 
approximation,  inspite  of  the  well-known  limitations  of  the  theory  which 
results  therefrom.  The  justification  for  this  severe  restriction  lies  in  the 
usefulness  of  simple  Hflckel  theory  in  correlating  diverse  properties  of  con¬ 
jugated  systems,  e.g.  Coulson  ,  especially  in  finding  out  the  biochemical 
potentiality  of  complex  molecules  and  as  a  starting  point  for  more 
sophisticated  theories  of  the  electronic  structure  of  these  molecules. 


C.A.  Coulson,  Symposium  on  Molecular  Physics,  Nikko,  Japan  (1953). 

It  is  of  great  interest  in  this  connection  to  note  that  recently  Lbwdin 
generalizing  certain  ideas  of  Brueckner  has  shown,  from  the  view -point 
of  perturbation  theory,  how  it  is  possible  to  obtain  from  the  eigenfunctions 
of  a  "model  Hamiltonian"  of  the  Hflckel  type  the  exact  eigenfunctions  of  the 
system  including  all  many -particle  interactions  exactly  by  means  of  the 
reaction  operator.  Hence,  in  principle,  we  can  pass  from  the  eigenfunctions 
of  the  Hdckel  Hamiltonian  to  the  exact  eigenfunctions  of  the  system,  at 
least  formally.  So  a  simple  and  direct  way  of  obtaining  the  eigensolutions 
of  the  Hilckel  Hamiltonian  for  conjugated  systems  would  be  of  considerable 
value. 


2b) 


P.O.  Ldwdin,  Technical  Notes  Nos.  47  and  48  from  the  Uppsala  Quan¬ 
tum  Chemistry  Group  (I960);  J.  Math.  Phys,  •••  (1962). 
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A,  straightforward  solution  of  the  eigenvector  and  eigenvalues  of 
the  Hflckel  Hamiltonian  for  complicated  molecules  could  be  rather  time- 
-consuming,  if  the  calculation  has  to  be  repeated  for  a  range  of  values 
of  the  basic  parameters  which  enter  into  the  theory.  By  a  perturbation 
approach,  we  can  use  maximally  the  already  known  information  of  the 
basic  skeleton  reducing  the  numerical  magnitude  of  the  problem  consider¬ 
ably  and  yet,  using  the  method  of  L.6wdin  obtaining  the  exact  solutions 
which  .correspond  to  the  use  of  the  whole  perturbation  series. 

In  section  I,  we  give  a  brief  review  of  some  of  the  important 
theoretical  work  done  on  the  theory  of  simple  conjugated  systems  in  the 
past.  Particular  attention  is  devoted  to  the  role  played  by  the  Hdckel 
theory  in  the  development  of  the  present  day  approach  to  the  calculations 
on  these  systems.  This  discussion  is  necessarily  limited,  and  it  is  meant 
only  to  set  the  stage  for  the  latter  sections.  Although  we  point  out  some 
of  the  successes  and  failures  of  the  simple  Hdckel  approach,  the  discussion 
is  definitely  not  a  critique  of  the  theory. 

In  section  II,  we  explain  the  two  methods  to  be  used  and  derive  the 
important  equations  which  we  shall  need  later  in  the  succeeding  sections. 

In  section  III,  we  illustrate  the  utility  of  the  first  or  non-iterative 
method  of  the  preceding  section  by  applying  it  in  the  calculations  of  the 
eigenfunctions  and  eigenvectors  of  some  locally  perturbed  ring  type  con¬ 
jugated  systems. 

In  section  IV,  we  continue  with  the  applications  and  demonstrate 
the  second  or  iterative  method  by  using  it  in  calculations  on  simple  con¬ 
jugated  systems  of  both  chain  and  ring  type. 

In  section  V,  we  conclude  by  giving  a  brief  discussion  of  the  two 
methods  amd  some  of  their  possible  applications. 
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I.  REVIEW  OF  WORK  ON  CONJUGATED  SYSTEMS 

In  1931,  Hflckel  applied  the  molecular  orbital  theory  developed 
by  Lennard- Jones,  Hund  and  Mulliken  in  an  investigation  of  the  ir-elec- 
tron  system  of  the  benzene  molecule.  Hdckel's  approach  to  the  problem 
can  be  characterized  by  two  basic  properties,  i.e.  the  well-known  (T-tt 
separation  of  the  electrons  and  the  two -parameter  nature  of  the  theory 
which  derives  from  the  use  of  the  so-called  "coulomb"  and  "i-esonance" 
integrals  usually  denoted  by  a  and  p  respectively.  Since  Hdckel's 
original  work,  a  substantial  body  of  literature  has  developed  around  the 
application  of  the  theory  to  conjugated  systems.  Hflckel's  theory,  especially 
in  its  semi -empirical  form,  has  been  an  extraordinarily  useful  tool  for 
the  correlation  of  mainly  ground-state  properties  of  conjugated  and  aromatic 
hydrocarbon  molecules.  The  theory  has  led  to  a  fundamental  understanding 
of  such  properties  as  resonance  stability,  bond  lengths,  charge  densities, 
free  valence  and  the  reactivity  of  these  molecules  to  electrophilic,  nucleo¬ 
philic  and  free  radical  reagents.  It  is  still  not  clearly  understood  why  the 
Hdckel  theory  is  so  successful  although  Hall  has  pointed  out  that  it  is 
most  likely  due  to  a  remarkable  correspondence  to  a  more  sophisticated 
approach  rather  than  to  the  validity  of  the  many  assumptions  inherent  in 
the  theory.  Ldwdin^^^  has  also  pointed  out  that  in  the  Hflckel  theory,  the 
Hamiltonian  operator  is  not  actually  specified  and  this  remarkable  success 
might  mean  that  the  two-parameter  Hamiltonian  matrix  may  correspond  to 

the  use  of  the  exact  effective  Hamiltonian  operator  including  all  many -body 

4b) 

interactions  exactly.  It  is  also  of  interest  to  mention  that  several  authors  ' 
have  pointed  out  that  indices  like  Coulson  charges  and  bond  orders  and 
polarizabilities  which  are  used  to  characterize  the  properties  are  precisely 
those  which  would  be  the  relevant  indices  for  the  molecule  from  a  considera¬ 
tion  of  the  one-particle  density  matrix  of  the  molecule.  So  certain  measure 
of  the  success  of  Hdckel  theory  is  due  to  the  fact  that  we  are  dealing  with 


E.  Hdckel,  Z.  PhysikTO,  204  (1931). 

G.G.  Hall,  Proc.  Roy.  Soc.  (London)  A213,  102,  113  (1952). 

P. O.  Lttwdin,  Proc.  Int.  Conf.  Theor.  Physics,  Japan,  Kyoto  1953, 
p.  599;  Phys.  Rev.  97,  1490  (1955);  R.  McWeeny,  Proc.  Roy,  Soc. 
(London)  A232,  114  (1955). 
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5) 

the  right  type  of  indices  for  these  molecules.  Pople's  calculations  '  on  the 
ground  state  of  aromatic  hydrocarbons  by  SCF  MO  theory  seems  to  justify 
this  conclusion.  Pople  made  a  number  of  approximations  which  had  the 
effect  of  making  the  SCF  LCAO  MO  equations  of  Roothaan  simpler  and 
more  amenable  to  solution.  He  also  showed  that  these  modified  equations, 
when  applied  to  alternant  hydrocarbons,  had  the  same  form  as  the  equations 
of  Hdckel  theory  plus  some  perturbation  terms  due  to  explicit  inclusion  of 
electron  interaction  and  that  the  SCF  charge  distribution  in  the  ground  state 
of  alternant  hydrocarbons  is  as  uniform  as  in  Hiickel  theory. 

The  applicability  of  the  Hiickel  method  can  be  extended  to  hetero- 

7) 

systems  '  like  pyridine  by  introducing  a  "coulomb  parameter"  appropriate 
for  the  substituted  atom  and  an  "inductive  parameter"  which  is  meant  to 
account  for  the  change  of  electronegativity  of  the  carbon  atoms  adjacent  to 

8) 

heteroatoms.  From  detailed  calculations  on  pyridine,  Wheland  and  Pauling  ' 

concluded  that  the  value  of  the  inductive  parameter  must  be  approximately 

one  tenth  of  the  coulomb  parameter  to  explain  the  orientation  in  reactions 

of  pyridine  and  in  later  works,  this  ratio  or  a  value  close  to  it  has  been 

used  extensively.  In  recent  years  and  in  light  of  more  refined  calculations 

9) 

of  an  SCF  type  such  as  those  reported  by  McWeeny  and  Peacock  '  and  by 
Brown  and  Heffernan  considerable  doubt  has  been  cast  on  the  validity  of 
using  the  inductive  parameter  originally  introduced  by  Wheland  and  Pauling. 

Coulson  and  Rushbrooke  ^  showed  that  the  Hdckel  theory  within 
its  initial  approximations  is  internally  self-consistent  when  applied  to 
alternant  hydrocarbons.  By  self-consistent  it  is  meant  that  the  assumption 
that  all  of  the  diagonal  elements  of  the  Hamiltonian  matrix  are  equal  leads  to 
the  conclusion  that  in  the  ground  state,  the  Coulson  charges  of  the  individual 
atoms  are  all  equal  and  thus  does  not  contradict  the  original  assumption. 


5) 

6) 

7) 

8) 

9) 

10) 
11) 


J.A.  Pople,  Trans.  Faraday  Soc.  1375  (1953). 

C.C.J.  Roothaan,  Revs.  Modern  Phys.  23,  69  (1951). 

E.  Hiickel,  Z.  Physik  ^2.  310  (1931);  628  (1932). 

G.W,  Wheland  and  L.  Pauling,  J.  Am.  Chem*  Soc.  57,  2086  (1935). 

R,  McWeeny  and  T.E.  Peacock,  Proc,  Phys.  Soc.  ALXX,  41  (1957). 

R.D.  Brown  and  M.L.  Heffernan,  Australian  J,  Chem.  10,  2ll  (1957). 

C.A.  Coulson  and  G.S.  Rushbrooke,  Proc.  Cambridge  Phil.  Soc.  36, 
193  (1910). 
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12a) 

Coulson  and  Longuet-Higlins  '  in  a  series  of  papers  generalized 
the  ideas  of  Wheland  and  Pauling  and  developed  a  second  order  perturba¬ 
tion  theory  for  conjugated  systems.  They  derived  some  very  basic  theorems 
concerning  the  total  charge  orders  and  bond  orders  of  the  xr-electron  system 
and  defined  certain  polarizability  coefficients.  These  coefficients  have 
proved  extremely  useful  in  understanding  the  chemical  reactivity  and  the 
changes  produced  in  the  properties  of  a  conjugated  molecule  on  substitution 
of  a  carbon  atom  by  a  heteroatom.  Lionguet-Higgins  has  applied  these 
polarizability  coefficients  in  an  understanding  of  the  basic  properties  of 
heteromolecules. 

13) 

Ldwdin  '  has  extended  the  theory  and  generalized  the  definition 
of  charge  orders  and  bond  orders  to  account  for  the  non -orthogonality  of  the 
basic  atomic  orbitals.  He  also  emphasized  the  importance  of  the  charge 
orders  and  bond  orders  of  the  individual  molecular  orbitals. 

The  choice  of  the  value  to  assign  for  the  coulomb  parameter  of 
the  substituted  atom  has  always  been  a  problem.  Wheland  and  Pauling  based 
their  choice  on  consideration  of  the  electronegativity  differences  between 
the  heteroatom  and  carbon  atom  and  both  the  Pauling  and  Mulliken 
scales  of  electronegativity  have  been  used.  It  is  evident  now  that  the  value 
of  the  coulomb  parameter  which  should  be  assigned  depends  on  the  physical 
quantity  under  consideration.  Thus  one  value  of  the  coulomb  parameter  will 
give  good  agreement  with  reactivity  data  whereas  a  rather  different  value 
is  necessary  if  agreement  with  dipole  moment  data  is  sought.  Illustrations 
of  this  fact  are  numerous  in  the  literature  but  we  might  mention  the  cal¬ 
culations  of  Longuet -Higgins  and  Coulson  Orgel  et  al  Lttwdin 


12a) 


12b) 

13) 

14) 

15) 

16) 

17) 

18) 


C.A,  Coulson  and  H.  C.  Longuet -Higgins,  Proc.  Roy.  Soc.  (London) 
A191.  39  (1947);  A192,  16  (1947);  A193.  447,  456  (1948);  A195,  188 
(1948). 

H.C.  Longuet-Higgins,  J.  Chem.  Phys.  18,  275  (1950). 

P.O.  Ldwdin,  J.  Chem.  Phys,  21,  496  (1953). 

L.  Pauling,  J.  Am.  Chem.  Soc.  3570  (1932). 

R.S.  Mulliken,  J.  Chem.  Phys.  2,  782  (1934);  3,  573  (1935). 

H.C.  Longuet-Higgins  and  C.A.  Coulson,  Trans.  Faraday  Soc. 

87  (1947). 

Orgel,  Cottrell,  Dick  and  Sutton,  Trans.  Faraday  Soc.  £7,  113(1951), 
Reference  13. 
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Lofthus  and  Hameka  and  Liquori  It  is,  however,  desirable 
theoretically  that  with  one  particular  value  of  the  parameter  we  should  be 
able  to  interpret  the  diverse  properties  of  heteromolecules. 

The  moat  recent  tendency  has  been  to  carry  out  some  type  of  SCF 

Cl  calculations  ibr  the  ir-electrons  with  integral  approximations  of  the 

2  ll 

type  introduced  by  Pariser  and  Parr  '  and  then  to  choose  the  parameters 
in  Hdckel  theory  in  such  a  way  that  the  SCF  results  for  the  charge 
densities  in  the  ground  state  are  reproduced.  The  calculations  reported 
in  references  9  and  10  are  of  this  type  and  Amos  and  Hall  have  reported 
a  similar  type  of  calculation.  The  results  have  led  to  the  conclusion  that 
the  value  of  the  coulomb  parameter  for  nitrogen,  for  example,  in  pyridine - 
-like  systems  should  be  between  0.2  and  0,4  in  order  to  reproduce  the 
SCF  results. 

The  simple  Httckel  theory,  so  useful  for  the  correlation  of  ground 
state  properties,  cannot  adequately  describe  excited  states.  Thus,  it  is 
not  possible  to  explain  the  spectra  of  benzene,  pyridine  and  the  disubstituted 
benzene  molecules  on  the  basis  of  the  naive  idea  of  transitions  between  the 
filled  and  unfilled  Hdckel  molecular  orbitals.  In  order  to  obtain  agreement 
with  the  experimental  results  on  the  spectra  of  these  molecules,  it  is 

Z3) 

necessary  to  carry  out  rather  elaborate  calculations.  Longuet -Higgins  ' 
has  discussed  the  excited  state  calculations  in  a  review  article  and  here 
we  shall  mention  only  two  of  the  problems. 

24) 

Coulson  '  showed  that  the  calculated  electronic  spectra  of  the 
polyacenes  could  be  made  more  consistent  with  the  calculated  ground  state 
properties  if  one  took  into  account  the  overlap  integral  between  neigh¬ 
bouring  carbon  atoms  instead  of  neglecting  it  as  Hdckel  originally  did. 


19) 

20) 
21) 
22) 

23) 

24) 


A,  Lofthus,  Molecular  Physics  9  (1958). 

H.F.  Hameka  and  A. M,  Liquori,  Molecular  Physics  367  (1959)« 
R.  Pariser  and  R.G,  Parr,  J»  Chem.  Phys,  21,  466,  767  (1953). 
A.T,  Amos  and  G."G.  Hall,  Molecular  Physics  4,  25  (1961). 

H.C.  Longuet-Higgins,  Advances  in  Chemical  Physics  1,  239  (1958). 
C.A.  Coulson,  Proc.  Phys,  Soc.  (London)  A68,  81  (1955). 
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The  inclusion  of  overlap  had  the  effect  of  distributing  the  energy 
of  the  unoccupied  molecular  orbitals  wider  and  compressing  the  occupied 
molecular  orbitals  in  energy  scale  which  leads  to  a  slightly  better  under¬ 
standing  of  the  spectra  of  polyacenes,  the  spectra  being  interpreted  to  be 
due  to  one  electron  jumps  from  occupied  Hd'ckel  molecular  orbitals  to 
unoccupied  onest  A  »tudy  of  the  problem  of  inclusion  of  overlap  was  made 
by  Ldwdin  and  later  by  Chirgwin  and  Coulson  and  the  former 

showed  that  overlap  could  be  included  by  carrying  out  a  symmetric 
orthogonalization  of  the  basis  functions  which  induces  a  linear  transforma¬ 
tion  of  the  Hamiltonian  matrix.  He  also  pointed  out  that  the  off-diagonal 
elements  of  the  transformed  Hamiltonian  matrix  are  invariant  to  shift  of 
the  origin  of  energy  while  those  of  the  original  overlap  neglected  Hamil¬ 
tonian  matrix  are  not,  and  therefore  any  parametrisation  of  the  Hamil¬ 
tonian  matrix  must  be  on  the  transformed  Hamiltonian  matrix,  rather  than 
on  the  original  Hamiltonian  matrix,  as  is  usually  done.  The.  case  of  hetero- 
atoms  has  further  been  studied  by  Del  Re  , 

Although  the  inclusion  of  overlap  can  improve  the  excited  singlet 
states,  the  problem  of  triplet  states  remains  until  one  attempts  to  take 
account  of  electron  repulsion  in  a  satisfactory  manner.  Unless  electron 
repulsion  is  included  in  the  calculations,  the  effects  of  overlap  cannot  be 
definitely  determined  as  Pople  has  shown. 

27) 

Goeppert-Mayer  and  Sklar  very  early  (1938),  tried  to  calculate 
the  spectra  of  benzene  using  a  Hamiltonian  which  included  repulsion  between 
the  ir-electrons.  Their  calculations  led  them  to  conclude,  correctly,  that 
the  excited  singlet  levels  were  separated  from  each  other  and  from  the 
triplet  levels  by  the  electron  repulsion  terms, but  numerical  errors  and 
the  necessity  of  approximating  many  center  integrals  prevented  them  from 
obtaining  accurate  quantitative  results.  Following  Goeppert-Mayer  and 
Sklar,  many  people  have  worked  on  the  theoretical  calculations  of  the 


P.O.  Ldwdin,  Arkiv  Math. ,  Fys, ,  Astr.  35A,  No,  9  (1947); 

J.  Chem.  Phys.  365  (1950);  Adv.  Phys.  5,  1  (1956). 

B.H.  Chirgwin  and  C.  A.  Coulson,  Proc.  Roy.  Soc.  (London)  A201, 
196  (1950). 

G.  Del  Re,  Technical  Note  No.  20  from  Uppsala  Quantum  Chemistry 
Group  (1958).  Published  in  Nuovo  Cimento  JT,  644  (I960). 

J.A.  Pople,  Trans.  Faraday  Soc.  1375  (1953). 

M.  Goeppert-Mayer  and  A.  L.  Sklar,  J.  Chem.  Phys.  6,  645  (1938). 


25a) 

25b) 

25c) 

26) 

27) 
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spectra  of  conjugated  molecules.  It  is  not  possible  for  us  to  give  an  account 

of  the  development  here  but  some  of  the  most  important  work  is  contained 

in  papers  by  Roothaan  and  Parr  Nira  Parr  et  al  Dewar  and 
31)  32) 

Longuet -Higgins  '  and  Pariser  Most  of  this  work  has  been  concerned 

with  the  problem  of  taking  into  account  electron  repulsion  (between  the 
individual  ir-electrons  and  between  the  a  and  ir  electrons)  in  a  satisfactory 
way. 


28) 

29) 

30) 

31) 

32) 


C.C.J.  Roothaan  and  R.G.  Parr,  J.  Chem.  Phys.  17,  1001  (1949). 

K.  Nira,  J.  Chem.  Phys.  1498  (1952). 

R.G.  Parr,  D.P.  Craig  and  I. G.  Ross,  J.  Chem.  Phys.  18,  1561  (1950). 

M.J.S.  Dewar  and  H.C.  Longuet-Higgins,  Proc.  Phys.  Soc.  (London) 
A^,  795  (1954). 

R.  Pariser,  J.  Chem.  Phys.  24,  250  (1956). 


We  quote  from  the  review  article  by  Longuet-Higgins  mentioned 

above: 


"It  may  be  said,  therefore,  that  the  electronic  spectra  of  aromatic 
hydrocarbons  are  now  well  understood.  The  essential  features  of  the  present 
theory  are 

(1)  that,  if  suitable  molecular  orbitals  are  adopted  in  the  determinant 
representing  the  wave  function  of  the  ground  state,  then  the  excited  states  may 
be  adequately  represented  in  terms  of  singly  excited  configurations  alone, 
provided  that, 

(2)  if  the  molecular  orbitals  are  not  self-consistent,  then  singly 
excited  configurations  contribute  to  the  ground  state  and  the  unexcited  conr- 
figuration  contribute  to  the  excited  states  and  provided  also  that, 

(3)  atomic  orbital  overlap  is  neglected  and  the  integrals  determining 
the  matrix  elements  of  the  Hamiltonian  are  evaluated  semi-empirically 
rather  than  by  nonempirical  calculation.  " 

We  have  not  yet  mentioned  either  the  method  of  alternant  molecular 

33) 

orbitals  (AMO)  suggested  by  Lowdin  '  or  the  method  of  split  p-orbitals 


33) 


P.O.  Ldwdin,  Nikko  Symp.  Mol.  Phys.  (Maruzen,  Tokyo  1954)  p.  13j 
Phys.  Rev.  97.  1509  (1955). 
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(SPO)  introduced  by  Dewar  and  Wulfman  in  the  former  method,  which 

is  applicable  to  systems  with  an  alterncUit  character,  electrons  with  plus 

spin  are  assigned  to  orbitals  which  have  large  amplitude  at  every  alternate 

atom  and  small  amplitude  at  the  remaining  atoms  while  electrons  with 

0 

minus  spin  are  assigned  to  orbitals  whose  amplitude  changes  alternately 
but  reversed  to  the  former  ones.  Thus  an  attempt  is  made  to  incorporate 
electron  correlation  directly  in  the  wave  function.  The  principal  idea  of 
the  SPO  method  is  quite  similar  in  that  electron  correlation  is  accounted 
for  by  constraining  the^  electrons  of  different  spin  to  be  in  different  parts 
of  space,  in  this  case,  the  two  lobes  of  the  p-type  atomic  orbitals.  It  may 
be  remarked  that  the  SPO  method  has  been  criticized  severely  by  Griffiths 

35) 

Itoh  and  Yoshizumi  '  have  applied  the  AMO  method  in  the  calcula¬ 
tion  of  the  grovind  state  energy  of  benzene  and  have  found  that  they  could 

account  for  85  o/o  of  the  correlation  error  as  determined  by  the  Parr,  Craig 
30) 

and  Ross  '  calculation  which  was  based  on  very  extensive  configuration 
interaction.  Recently,  Pauncz,  de  Heer  and  Lttwdin  derived  a  general 
energy  expression  for  an  arbitrary  alternant  system  with  closed  shell 
structure  and  made  a  numerical  analysis  of  cyclic  systems  consisting  of 
2n  =  4v  +  2  electrons  moving  in  the  field  of  2n  C-H  centers.  They  arrived 
at  the  result  that  the  energy  depression  decreases  only  slowly  with  increas¬ 
ing  n  and  hence  the  method  should  be  useful  for  large  systems.  Thus,  it 
appears  that  the  AMO  method  may  indeed  be  a  very  simple  device  for  in¬ 
cluding  electron  correlation  in  a  particular  class  of  molecules.  The  SPO 
method  has  been  applied  to  ethylene  by  Dewar  and  Hojvat  with  excellent 
results  and  calculations  are  reported  to  be  underway  on  the  ground 
state  of  benzene.  It  will  be  interesting  to  compare  the  results  of  the  AMO 
and  the  SPO  calculation  on  this  molecule. 


34a) 

34b) 

35) 

36) 


37) 


M^J.S.  Dewa,r  knd  C.E*  Wulfman,  J.  Chem.  Phys.  158  (1958). 

J.S.  Griffiths,  J.  Chem.  Phys.  36,  1689  (1962). 

T»  Itoh  and  H.  Yoshizumi,  J,  Phys.  Soc.  Japan  10,  201  (1955). 

R.  Pauncz,  J.  de  Heer  and  P. O.  Lbwdin,  Technical  Notes  Nos.  55 
and  56  from  Uppsala  Quantum  Chemistry  Group  (I960);  published 
in  J.  Chem.  Phys.  2247,  2257  (1962). 

M.J.  S.  Dewar  and  NX.  Hojvat,  J.  Chem.  Phys.  34,  1232  (1961). 
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Actually,  it  seems  to  be  quite  possible  that  the  two  methods  can  be 
and  should  be  combined  (CAMO  SPO)  for  the  following  reason.  The  AMO 
method  accounts  for  correlation  in  the  direction  around  the  ring  whereas 
the  SPO  method  accounts  for  the  correlation  in  a  direction  perpendicular 
to  the  plane  of  the  molecule.  Thus  the  situation  here  is  analogous  to  the 
radial  and  angular  correlation  in  atoms.  Both  types  of  correlation  should 
be  included  simultaneously  in  any  calculation  if  one  wishes  to  determine 
the  importance  of  one  type  compared  to  the  other.  It  would  seem  to  the 
authors  that  the  calculations  on  the  benzene  molecule  with  the  AMO  approach 
probably  over -emphasize  the  around -the -ring  correlation. 

Both  the  AMO  and  the  SPO  methods  can  be  considered  as  special 
cases  of  the  general  idea  of  different  orbitals  for  different  spins  which  has 
been  discussed  frequently  in  the  literature 


P.O.  Lbwdin,  Phys.  Rev.  97^,  1509  (1955);  Adv.  Phys.  1  (1956). 


We  shall  conclude  this  section  by  giving  a  summary  of  the  physical 
properties  of  interest  in  conjugated  systems  and  the  success  or  lack  of 
success  which  simple  Hdckel  theory  has  had  in  dealing  with  them. 

a)  Bond  lengths. 

The  bond  length  is  correlated  with  the  calculated  total  ir-bond  order, 
p  ,  which  is  defined  by 

X  o 

tr-S  "  ^  (1) 

Here  j  runs  over  the  occupied  molecular  orbitals,  m  is  the  occupation 
number  of  the  j-th  orbital  and  the  c's  are  the  coefficients  of  the  oribtals 
in  the  basis  set.  If  r  and  s  are  not  neighbours  the  bond  order  is  taken  to 
be  zero.  The  correlation  is  carried  out  by  drawing  a  bond  length  vs  bond 
order  curve  based  on  several  experimentally  known  points  for  a  particular 
type  of  bond.  This  curve  is  then  used  to  predict  the  bond  length  in  other 
molecules  containing  the  same  type  of  bond.  The  Hdckel  theory  is  really 
quite  successful  in  predicting  the  bond  lengths  of  conjugated  emd  aromatic 
molecules  to  an  accuracy  of  about  0.02  A  which  is  sufficient  for  most 
practical  purposes  besides  being  almost  at  the  limit  of  the  accuracy  of  two 
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dimensional  X-ray  analysis. 


Detailed  calculations  by  Coulson  et  al 


39) 


C.A.  Coulson,  R.  Daude!  and  J.M.  Robertson,  Proc.  Roy.  Soc. 
A207,  306  (1951). 


have  shown  that  configuration  interaction  and  other  refinements  do  not 
greatly  alter  the  relative  value  of  bond  orders.  The  bond  order  is  not  very 
sensitive  to  the  form  of  the  wave  function  and  thus,  inspite  of  its  great 
value  in  other  respects,  cannot  be  used  as  a  criterion  for  a  satisfactory 
wave  function. 

b)  Charge  Densities. 

The  ir-electron  charge  density  on  any  atom  in  a  conjugated  molecule 
is  correlated  with  the  total  charge  order  which  is  given  by 

%.=■  P) 

where  the  symbols  have  the  same  meaning  here  as  in  equation  (1).  The 
charge  orders  and  bond  orders  form  the  elements  of  the  charge  and  bond 
order  matrix  which  for  these  systems,  is  twice  the  first  order  density 
matrix  defined  by 

P=  (3) 

The  charge  orders  calculated  with  Hfickel  theory  are  sufficiently  accurate 
to  enable  the  chemist  to  predict  the  most  probable  position  at  which  a  charge 
reagent  will  attack  the  molecule. 


P.O.  Lbwdin,  Phys.  Rev.  97^,  1474  (1955). 

R.  McWeeny,  Revs.  Modern  Phys.  32,  335  (I960). 
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c)  Electric  dipole  moment. 

The  ir-charge  distribution  calculated  by  Hdckel  theory  can  yield 
a  qualitative  understanding  of  the  values  of  dipole  moments  among  a 
series  of  molecules.  In  order  to  obtain  correlation  with  experimental  data 
it  is  necessary  to  choose  values  for  the  parameters  which  differ  consider¬ 
ably  from  those  values  which  give  the  best  results  for  other  physical 
properties.  One  must  be  careful,  when  treating  dipole  moment  data,  to 
take  account  of  the  contributions  from  lone  pair  electrons  and  from  the 
O’ -framework.  The  fact  that  the  dipole  moment  is  quite  sensitive  to  small 
changes  in  the  wave  function  which  leave  the  energy  almost  unaffected 
further  complicates  the  problem. 

d)  Orbital  energies  and  spectra. 

The  simple  Hdckel  theory  is  incapable  of  giving  even  satisfactory 
qualitative  agreement  with  the  experimental  data  on  the  spectra  of  con¬ 
jugated  systems. 
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II.  DERIVATION  OF  EQUATIONS 

We  consider  a  system  whose  Hamiltonian  operator  can  be  written  as 

H  —  Ho  4-\/  (4) 

and  whose  exact  eigenfunction  ¥  can  be  expanded  in  terms  of  a  set  of  basis 
functions  basis  set  can  be  chosen  as  orthonormal  i.e. 

=  Sfet,  W 

We  can  carry  out  the  usual  transformation  of  the  original  eigenvalue  equation 
H?  =  EY  into  the  matrix  equation 

((HLtV)c 

in  which 

and  C  is  the  column  matrix  whose  elements  are  the  expansion  coefficients 
of  Y  in  the  basis  set  Equation  (6)  can  be  rewritten  in  the  form 

C=  (ED (8) 

Following  L6wdin  in  reference  I,  we  define 

O  ~  "  W*)  (9) 

so  that 

C?=G\/^(C  (10) 

oC,^ 

If  S’  can  be  obtained  in  closed  form  and  if  is  a  matrix  with  only  a 

few  nonvanishing  elements,  some  rather  simple  expressions  for  E  can  be 
derived.  As  an  example,  we  can  consider  the  case  when  only  one  element  is 
non-vanishing,  say,  ®  • 


EC 


(6) 


<SfelV|S4> 


(7) 
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Equation  (10)  then  becomes 

C,  =^,V„et  or  G)i  '^11  '"f  (11) 

and 

~  (12) 

Equation  (11)  has  been  used  by  Koster  and  Slater  and  by  Lax  The 


42) 

43) 


G.F.  Koster  and  J.C.  Slater,  Phys.  Rev.  1167  (1954). 
M.  Lax,  Phys.  Rev.  94.  1391  (1954). 


former  authors  evaluated  the  matrix  by  means  of  the  transformation  of 

Ho  to  diagonal  form,  whereas  Lowdin  has  pointed  out  that  there  are  many 
alternative  ways  for  evaluating  this  inverse  matrix  and  that  actually  in  many 
problems,  the  matrix  itself  is  not  needed  but  only  the  vector  solution  of  a 
linear  equation  system. 

The  systems  which  we  shall  consider  here  all  have  the  characteristic 

property  that  within  the  framework  of  simple  Hflckel  theory,  the  matrix  (^0^ 

is  cyclic  or  can  be  chosen  as  cyclic  with  suitable  choice  of  the  perturbation 

matrix  .  It  is  customary  to  use  only  the  "tight -binding  approximation"  or 

"nearest  neighbour  interaction",  which  we  shall  also  use.  This  cyclic  nature 

of  allows  us  to  find  (fy  by  using  the  powerful  method  developed  by 

44) 

Lfiwdin,  Pauncz  and  de  Heer  '  in  a  closed  analytic  form  in  terms  of 
Chebyshev  polynomials 


44) 


45) 


P.O,  Lfiwdin,  R.  Pauncz,  and  J.  de  Heer,  J.  Math,  Phys.  1,  46l 
(I960):  It  is  interesting  to  note  that  E.  Bodewig  has  used  a  similar 
technique  (Matrix  Calculus  by  E.  Bodewig,  North  Holland  Publishing 
Co.,  Amsterdam).  However,  he  does  not  mention  the  fact  that  the 
polynomials  which  he  introduces  are  actually  nothing  but  the 
Chebyshev  polynomials. 

See  appendix  for  a  summary  of  the  properties  of  these  polynomials 
and  refer  to  the  "Tables  of  Chebyshev  polynomials"  published  by 
National  Bureau  of  Standards  in  their  applied  Mathematics  Series  for 
a  more  detailed  discussion. 
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We  shall  now  give  a  brief  description  of  this  naethod.  The  physical 
model  which  we  have  in  mind  during  the  discussion  is  that  of  a  number^  n^ 
of  atoms  placed  equidistant  from  each  other  on  the  circumference  of  a  circle* 
Interactions  between  nearest  neighbours  only  are  considered.  We  must 
distinguish  two  cases  a)  an  even  number  of  atoms  n  =  2v  and  b)  an  odd  num¬ 
ber  n  =  2v  +  1  . 

a)  n  =  2v 

In  this  case,  we  use  the  equation  (equation  34  of  reference  44) 

The  ,  which  are  cyclic  and  symmetric  will  be  called  topological 

matrices  after  Ruedenberg  who  has  discussed  I  .  The  matrix 
elements  of  ,  which  we  shall  relate  to  our  problem  below,  are  given 

by 


if  p  and  q  are  neighbours 
otherwise 


(14) 


It  is  somewhat  easier  to  manipulate  the  if  these  are  expressed  in 

terms  of  another  type  of  elementary  matrix.  We  write 

M-fc.  =  (») 

where 

yn^  =  w 

is  the  transpose  of 


K.  Ruedenberg,  J.  Chem.  Phys.  34,  1884  (1961). 


The  elements  of  the  yn  matrices  are  given  by 


46) 


The  prime  is  to  remind  us  of  the  cyclic  nature  of  the  problem  and  simply 
means  that  if  (p  +  k)  is  outside  of  the  range  (p  +  k)^  n  we  add  n  if 
(p  4-  k)  <  1  and  subtract  n  if  (p  4-  k)  >  n  ,  e.g>  if  n  =  10  ,  k  =  2  and 
p=9  then  6(n).,i  =  S,i‘ 

The  are  the  Chebyshev  polynomials  of  the  first  kind  (see 

appendix)  and  (a)  is  defined  by 

FvW  =•  (18) 

Also,  6  =  c  =4-  and  all  other  e,  =  1  . 

’  o  V  2  k 


b) 


n  =  2v  4-  1 


47) 


Here  we  use  the  equation  (equation  43  of  reference  44) 


(19) 


We  use  the  Chebyshev  polynomials  of  the  first  kind  in  this  formulation 
but  it  is  often  more  convenient  to  use  the  polynomials  of  the  second 
kind  as  we  shall  do  in  Section  IV. 


with 


Gv  “  ^ y-f/  ^ 


(20) 


3.11  other  qj^  (including  q^  )  are  equal  to  1. 

According  to  the  assumptions  of  simple  Htickel  theory,  fa)^4=  ■ 

f  element,  of  ate  teto. 

a  and  p  are  parameters  which  are  usually  determined  empirically.  U  sing 
these  quantities  we  can  write 
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where 


Equation  (8)  now  becomes 


e=  —h  (n,+a-i)  VC 


This  equation  can  be  made  much  simpler  if  the  energy  is  expressed  in  units 
of  p  •  Then  equation  (23)  becomes 


c  ~  'v(C 


(23a) 


Now  from  equation  13,  we  have 


This  can  be  written  in  terms  of  the  T^T^^niatrices,  so  that 
Using  the  properties  of  the  "l^)^J^matrice8,  this  becomes 

C&vlci  =  t  ]  (26) 

'hmo  V-  Y—l  ■'  ^ 


Equation  (10)  now  becomes 


(■  ^ (W,0  Hi  e,  (2’) 

4/  ^->2  1—1 


and  finally 

This  is  the  general  equation  for  the  case  n  =  2v  which  we  shall  use. later. 
The  corresponding  equation  for  the  case  n  =  2v  +  1  is 

Ci  =  -  i  %  X-*  w  n  %+^oy  +  V  <y  (2’) 
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These  last  two  expressions  give  rise  in  application  to  equation  systems 
whose  order  depends  upon  the  number  of  nonvanishing  matrix  elements  in 
the  perturbation  matrix  \/J  .  The  equation  system  is  easily  solved  for  one 
and  two  non -vanishing  elements  and  frequently  the  symmetry  of  the  problem 
makes  it  possible  to  solve  a  set  of  equations  of  order  three.  More  compli¬ 
cated  perturbation  matrices  give  rise  to  equation  systems  which  are 
treated  better  by  a  second  order  iteration -variation  procedure  We 

shall  now  give  a  brief  description  of  this  method. 


48) 


P.O.  Lfiwdin,  Technical  Note  No.  11,  Uppsala  Quantum  Chemistry 
Group  (1958),  and  the  reference  to  the  previous  literature  given 
there. 


Following  Lfiwdin  we  start  from  equation  10  and  rewrite  it  as 


(l-(&v)(r  =  0 


(30) 


We  note  that  (■D  -  G-V )  is  non-hermitiam.  Since  it  is  easier  to  apply 
partitioning  techniques  to  a  hermitian  matrix,  we  multiply  by  V  from  the 
left  and  obtain 


(V  -  VGV)a:  =  0 

where  (Ns 


0Y»  IN1(D  =  0 


(31) 

(32) 


Next  we  partition  the  basis  set  (  4>,  }  into  two  subsets  4i  ,  <l)o  . . .  $ 

and  from  equation  (31)  obtain  the  two  equations  I  ? 


Nn  C|  +  —  0 

(33) 

N41 C,  +•  =  0 

(34) 

From  these  two  equations  we  then  obtain 

(N|,  -  Nn  m;;  c,  ^0 

(35) 

If  Cj  0 


N„  -  IWk,  M)^^'  =  F(e)  r:  o 


(36) 


This  equation  contains  E  implicitly  and  is  best  solved  iteratively.  The 
iteration  procedure  as  shown  by  Ldwdin  in  reference  48  corresponds  to 
the  application  of  the  Newton -Raphson  process  to  F(E)  .  Thus 


where  the  prime  denotes  differentiation  with  respect  to  E  .  Now 


(37) 


F'(e=  N,1  -  WuHtK.  ■  (3“) 


We  can  simplify  this  equation  considerably  by  using  equation  (34)  and  its 
complex  conjugate  so  that 


r'teo  =  N,;  +  |N),(C1  +  ®t,  w/j  Ct  ^ 

c,  -^js- 

which  can  be  written,  very  simply,  as 

P^E}  =  €*H'(C 

|C.|^ 

Also 

F(e1—  /^u  -f  IMi/.<r-4  —  ('Cm  Cl  +  lv)|4(r{  „ 

c,  57-^  - 

where  the  subscript  means  we  take  the  first  element  of  jpcj 
fore,  we  have  that 

[:w(^,  ie,r 


E*=  E 


Now  from  the  definition  of  ^  and  we  obtain 


(39) 


(40) 


(41) 

.  There- 

(42) 


W  =  V  (  EtJ  -  IH).)1 7} .  (  S'!)  -  M.)  V 

“  V  C"  C-  V  —  V 

and  therefore, 

t  =•  b  ^ - - — 


(43) 


(44) 
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So  the  iteration  procedure  consists  choosing  a  trial  value  of  E  •  say 
E  ,  (different  from  the  eigenvalues  of  (H)e  )  and  solving  the  equation 


system  (from  equation  34) 


e«  = 


(45) 


If  we  '■hoose  c^  =  1  , 


then 


c  = 


to 


(46) 


We  get 


E.  - 


Cv('ll-G.V)4g 


where 


<&o=  (E.'I-iH}.) 

define  (Gv  C^e 


(47) 


(48) 


We  define  (Lr*  V'Cen.  CLo  and  notice  that  for  the  exact  E  , 

unite  the  second  order  iteration  formula  as 

'  „„ 

<Iao  Clo 

In  the  perturbation  of  homocyclic  conjugated  systems  to  heterosystems, 
the  perturbation  usually  involves  only  a  few  elements  in  the  matrix, 
so  that  the  equation  system  to  be  solved  will  be  of  smaller  order  than  that 
of  the  original  problem.  For  example,  for  perturbations  on  three  atoms, 
the  equation  system  Nc  =i0  would  give  only  three  equations  involving 
three  coefficients  only,  and  one  of  these  coefficients  can  be  determined  from 
the  normalization  condition  on  •  This  whole  process  is  convenient  if  the 

matrix  (&  can  be  found  easily.  Fortunately,  all  heterosystems  both  open 
chain  and  cyclic,  can  be  achieved  by  suitable  perturbations,  diagonal  or  off- 
-diagonal  or  both  on  a  cyclic  system.  As  we  have  seen  above,  the  powerful 
method  of  L&wdin,  Pauncz  and  de  Heer  is  available  for  obtaining  for 

cyclic  systems  and  it  is  this  fact. that  leads  to  enormous  simplification  for  the 
perturbation  treatment. 
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m.  PERTURBED  RING  SYSTEMS 


In  this  section  we  shall  treat  various  locally  perturbed  ring  systems, 
particularly  the  mono-,  di-  and  tri- substituted  benzene  ring  by  the  first 
i. e.  the  non-iterative  method  outlined  in  the  previous  section.  Here  we 
shall  consider  only  simple  perturbations  with  two  and  three  non-zero 
elements  at  most  in  the  perturbation  matrix.  We  shall  neglect  the  non- 
-diagonal  inductive  effects  for  the  present  but  in  a  later  section  we  shall 
show  how  they  can  be  taken  into  account,  approximately  by  first  order  per¬ 
turbation  theory  within  the  framework  of  the  non-iterative  method.  It  is 
not  our  intention  to  relate  our  results  to  the  physical  properties  of  any  real 
molecule;  we  wish  only  to  illustrate  the  method  and  to  set  forth  certain 
relationships  and  data  which  we  believe  may  be  useful  in  later  studies.  We 
shall  treat  the  mono -substituted  six  membered  ring  rather  thoroughly  and 
only  indicate  how  some  of  the  results  can  be  obtained  for  the  more  complicated 
perturbations  and  systems. 


For  reference  purposes,  we  give  first  the  Hiickel  theory  result  for 
benzene.  There  are  two  non -degene rate  levels  given  by  a  +  2p  and  two 
doubly  degenerate  levels  given  by  a  +  P  where  a  and  P  are  the  coulomb 
and  exchange  parameters,  respectively.  The  lowest  three  levels  are  doubly 
occupied  and  the  other  levels  unoccupied  in  the  ground  state.  The  eigenvalues 
and  eigenvectors  are  shown  in  Table  I.  The  ^'s  are  the  atomic  orbitals  and 
the  numbering  scheme  is  given  by 


It  should  be  kept  in  mind  that  p  is  negative  and  therefore,  a  2p  is  the 
lowest  level. 


TABLE  L  Eigenvalues  and  eigenvectors  of  benzene. 


Level 

Energy 

Eigenfunction 

Symmetry 

C 

o  +  2p 

{6)'^(1  11111) 

^2u 

B 

o  +  P 

(12)'2(2  1  i  2  i  I) 

1 

®lg 

A 

a  +  P 

(4)’2 (0  1  1  0  I  i) 

1 

A' 

o  -  P 

{i2)"2{2  i  i  2  I  i) 

®2u 

B' 

a  -  p 

(4)“2  (0  1  I  0  1  I) 

®2u 

C 

CM 

1 

(6)'^  (1  I  1  i  1  i) 

_ 2g _ 
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T he  symbols  at  the  right  under  "Symmetry"  are  taken  from  the  usual 
group  theoretical  notation  whereas  the  symbols  to  the  left  under 
"Level"  represent  a  notation  in  frequent  use  in  the  literature  of  alter¬ 
nant  systems.  For  a  general  alternant  the  occupied  molecular  orbitals 
are  labelled  A,  B  ...  going  toward  larger  negative  energies  from  the 
energy  zero;  the  unoccupied  molecular  orbitals  A%  B'  ...  are  sym¬ 
metrically  located  on  the  positive  side  of  the  zero  of  energy  and  any 
pair  K  and  K'  differ  by  a  sign  change  in  the  coefficients  of  one  set 
of  atomic  orbitals.  The  same  notation  can  be  used  for  the  substituted 
alternant  if  we  agree  that  K  will  mean  that  orbital  which  goes  over 
into  K  of  the  parent  hydrocarbon,  as  the  substituents  regarded  as 
perturbation  are  removed.  This  convention  does  not  imply  that 
. . .  B,  A,  A',  B*  . . .  be  in  ascending  energy  order  as  they  are  in  the 
unsubstituted  case.  We  shall  use  this  notation  later  when  we  give 
diagrams  of  the  energy  levels  of  the  substituted  ring  and  we  shall  see 
that  two  levels,  as  a  function  of  the  perturbation,  can  actually  cross. 


a)  Mono -substituted  benzene  ring. 

If  we  consider  a  perturbation  matrix  V  with  a  single  non-zero 
element,  say  we  have  from  equation  (28), 

From  this  we  obtain,  on  carrying  out  the  summation 


^*Tv(*o  ^1,7]  o'l 


(51) 


or 


M,  = 


(52) 


It  is  possible,  of  course,  to  write  the  out  explicitly  so  that  wo 

have  an  implicit  equation  for  the  energy  parameter  a  in  terms  of  the 
perturbation  matrix  elements.  Thus  far,  we  have  found  it  more  practical 
and  instructive  in  the  cases  we  have  dealt  with,  to  use  tables  of 
Chebyshev  polynomials  directly  and  calculate  Vjj  for  a  number  of 
different  values  of  a  .  The  energy  E  can  then  be  graphically  displayed 
as  a  fvuiction  of  V  . 


Let  us  consider,  for  example,  a  six  membered  unperturbed  benzene 
ring.  The  substitution  of  a  nitrogen  atom  for  a  carbon  and  hydregon  atom 
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pair  gives  pyridine  as  our  perturbed  system  and  equation  (52)  becomes 


\/,=  _  = 
-CjW 


(53) 


It  is  a  matter  of  a  short  time  to  find  the  behaviour  of  all  of  the  energy 
levels  as  a  function  of  the  strength  of  the  perturbation.  We  show  this 
behaviour  in  Figure  1.  We  see  that  two  of  the  energy  levels  of  benzene 
remain  unperturbed  by  a  single  diagonal  perturbation.  This  is  to  be 
expected  since  there  are  two  molecular  orbitals  in  the  Hdckel  theory  of 
benzene  for  which  the  coefficients  of  the  atomic  orbital  centred  on  the 
site  of  the  localized  perturbation  vanish  and  overlap  and  inductive  effects 
have  been  neglected. 


The  asymptotes  for  the  energy  levels  can  be  found  simply  by  putting 
the  denominator  equal  to  zero  in  equation  (53).  that  is  to  say,  the  roots 
of  ~  ®  are  the  asymptotes  provided,  of  course,  that  they  are 

not  also  roots  of  the  numerator.  In  the  limit  of  very  large  perturbations, 
the  levels  which  split  off  are  given  by  V(a)  =  a  .  A  more  general  discus¬ 
sion  of  this  effect  has  been  given  by  Lbwdin  in  reference  1.  It  is  interest¬ 
ing  to  observe  that  if  j  =  0  we  have  Fj  =  0  which  gives  the  energy 
levels  of  the  unperturbed  problem.  Thus  we  can  always  write  F^  in 
factorized  form  since  we  will  always  know  the  unperturbed  energy  levels, 
e.g. 


=  (a-j-O  C^-O  [  CK+A.) 


(54) 


It  is  also  interesting  to  note  that  the  energy  levels  split  off  from 

the  top  and  bottom  of  the  "band"  just  as  in  the  case  of  localized  perturba- 
49) 

tions  in  solids  We  shall  see  later  that  we  have  one  level  splitting 

for  each  substituted  atom.  We  also  observe  that  the  energy  levels  which 
are  perturbed  all  move  in  the  same  direction. 

49) 

'  See  the  discussions  in  references  42  and  43. 


The  eigenvectors  of  the  perturbed  system  can  be  found  from  equation 
(28).  Thus,  we  find  for  C,^ 


The  energy  levels  of  pyridine  as  a  function  of  the  strength  of 
a  single  diagonal  perturbation.  Side  naarkers  indicate  the 
asymptotes.  A,  A'  etc.  label  the  symmetry  levels  as  dis¬ 
cussed  in  the  text  for  benzene. 
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But  from  equation  (53), 

so  that 

Cj.=  -'CiX^o. 

In  a  similar  way,  we  find  that 

c,  = 

C,  = 

^5"  “^3  Cin\Ji  ^stCa, 


(56) 

(57) 

(58) 

(59) 

(60) 


Of  course,  we  know  already  the  eigenvectors  going  with  the  unperturbed 
levels. 


To  find  the  eigenvectors  corresponding  to  any  value  of  the  perturba¬ 
tion  V7||  ,  it  is  only  necessary  to  use  Figure  1  to  pick  out  the  corre¬ 

sponding  eigenvalue,  turn  to  the  tables  of  Chebyshev  Polynomials,  tabulate 
the  If's  and  carry  out  the  operations  indicated.  The  calculation  of  eigen¬ 
vectors  for  any  value  of  Wn  can  be  carried  out  quickly,  for  this  case 
and  the  accuracy  is  limited  only  by  the  accuracy  with  which  the  eigen¬ 


value  can  be  read  from  the  graph  and  by  the  accuracy  with  which  the 
polynomials  are  tabulated  (3  decimal  places  in  the  tables  cited  in  reference 
45  with  an  interpolation  scheme  if  more  figures  are  needed).  For  future 
reference,  we  give  the  matrix  C  for  the  mono- substituted  benzene  ring 
when  =  1 


0.646 

0.517 

0.  000 

0.497 

0.000 

0.262 

0.413 

0.082 

0.500 

-0.423 

0.500 

-0.379 

0.295 

-0.409 

0.500 

-0. 198 

-0. 500 

0.455 

0.259 

-0.621 

0.000 

0.563 

0.000 

-0.481 

0.295 

-0.409 

-Oi 500 

-0. 198 

0.500 

0.455 

0.413 

0.082 

-0. 500 

-0.423 

-0.500 

-0.379 

Once  d?  is  known,  the  first  order  density  matrix  given  by  equation  (3) 
can  be  calculated  readily  and  since  we  also  may  wish  to  have  it  for 
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reference  in  later  work  we  give  it  here: 


<P 


0.685 

0.310 

-0.020 

-0.  154 

-0.020 

0.310 

0.310 

0.427 

0.338 

0.056 

-0. 162 

-0.072 

-0.020 

0.338 

0.504 

0.330 

0.004 

-0.  162 

-0. 154 

0.056 

0.330 

0.453 

0.330 

0.056 

-0.020 

-0. 162 

0.004 

0.330 

0.504 

0.338 

0.310 

-0.072 

0.  162 

0.056 

0.338 

0.427 

b)  Pi-subBtitnted  benzene  ring. 


Here,  we  consider  the  case  that  the  perturbation  matrix  V 
two  non-zero  elements  (assumed  equal),  both  located  on  the  diagonal.  From 
equation  (28),  two  relationships  among  the  can  be  obtained  and  the 

solution  of  these  gives  an  expression  for  the  strength  of  the  perturbation  in 
terms  of  Chebyshev  polynomials  of  the  energy  parameter.  Also  the  eigen¬ 
vectors  and  hence  the  density  matrix  can  be  readily  obtained  from  equation 
(28).  The  procedure  is  exactly  the  same  as  in  the  mono -substituted  case 
except  that  now  there  are  two  equations  instead  of  one  to  solve  so  that  the 
calculation  is  somewhat  more  involved  but  still  quite  simple. 


For  the  six  membered  ring  there  are  three  distinct  possibilities  for 
the  position  of  the  two  heteroatoms.  These  are  shown  in  the  following 
diagrams  with  the  corresponding  equation  for  V  and  the  eigenvectors 
(unnormalized) . 


PYRIDAZINE 


Ca  =  ±c,  ■  Cj=  •  Cqr  =  ±.^5- 

PYRIMIDINE 


V=  - 


c,+=:-^c<+  ) 

Cl  z%.  T» 


(63) 

(64) 


(65) 

(66) 


(67) 
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C| 


-  ±  Gj  =  ±  Cy 


4:  C4  ^  C<,  ±  G, 


(68) 

(69) 


The  simplicity  of  these  equations  is  such  that  it  is  possible  to  obtain  E  = 

=  E(V)  curves  in  a  short  time.  Care  must  be  taken,  however,  to  associate 
the  eigenvalues  and  eigenvectors  correctly.  In  the  above  equations,  the 
top  sign  in  the  denominator  of  the  energy  expression  always  must  be 
associated  with  the  top  sign  in  the  expression  for  the  coefficients.  For 
pyrimidine  we  have  introduced  an  obvious  notation  which  can  be  used  in 
lieu  of  the  above  convention  when  necessary.  It  may  seem  curious,  that, 
in  each  of  the  above  cases,  we  have  two  expressions  for  V  each  apparent¬ 
ly  of  the  fourth  degree  and  therefore  we  should  expect  eight  energy  levels 
altogether.  If  we  were  to  write  the  denominators  and  numerators  out  in 
terms  of  the  a's  we  would  find  that  in  all  cases  just  enough  roots  of  the 
unperturbed  problem  are  contained  in  the  denominator  to  cancel  the  corre¬ 
sponding  roots  in  the  factored  form  of  the  numerator  and  thus  reduce  the 
degree  of  the  equations.  For  example  we  have  for  pyridazine  in  terms  of  a 

Vice)  -  ^  Ca-»-0  (<^-0 

This  point  is  discussed  in  more  detail  and  for  the  general  case  in  reference  1. 

The  E(V)  vs  V  curves  for  the  disubstituted  six  membered  ring  are 
given  in  Figures  2-4.  We  note  that  two  levels  split  off  from  the  band  in 
these  cases.  The  eigenvectors  and  hence  the  density  matrices  can  be  ob¬ 
tained  in  a  straight  forward  manner. 


c)  Tri- substituted  benzene  ring. 

For  this  case  we  have  three  diagonal  elements  in  the  ^  matrix. 

As  in  the  case  of  the  di-substituted  ring,  there  are  three  distinct  possibilities 
for  the  positions  of  the  heteroatoms,  i.e. 


Figure  3,  The 


as  a  function  of  the  strength 
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We  9hall  consider  only  the  last  of  these  for  two  reasons:  1)  it  is  the  most 
stable  molecule  of  the  three  and  2)  the  symmetry  of  the  molecule  makes 
it  easy  to  solve  the  equations  by 'our  present  method.  We  find  that 

Here  we  have  two  expressions  with  different  terms  in  the  denominator 
rather  than  the  4-  signs  familiar  from  the  disubstituted  cases.  Also  the 
second  expression  here  is  the  same  as  the  expression  for  V  in  the  case 
of  pyrimidine  if  we  take  the  minus  sign  there.  This  is  easy  to  understand 
from  the  symmetry  of  the  molecule  and  the  form  of  the  unperturbed 
molecular  orbitals. 

For  the  eigenvectors  we  find  that  for  (a)  we  have 

0,=:C^=o  j  Cs=.-C5-  C  >  Ca  =  (72) 

T>- t,| 

For  (b)  we  have 

C5=.cr=c,  1  Ci  =  Ci  =  C4  pjj 

The  energy  curves  are  given  in  Figure  5  .We  note  that  three  levels 
split  off  in  this  case  although  two  of  them  remain  degenerate; 

d)  Perturbations  on  a  seven  membe red  ring. 

We  can  derive  equations  similar  to  the  above  for  seven  and  five 
membered  rings.  The  procedure  is  so  similar  that  there  is  no  point  in 
discussing  it  here.  For  example,  with  a  single  diagonal  perturbation 
equation  (29)  gives 

V'l,  -  ___ - — 

Dv  C«\) 

in  close  analogy  to  the  result  for  the  six  membered  ring. 


(74) 
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The  energy  levels  of  the  1,  3,  5  substituted  benzene  ring  as  a 
function  of  the  strength  of  the  perturbation. 
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®)  Perturbation  oji  a  ten  xnembered  ring. 

Both  Naphthalene  and  Azulene  can  be  considered  as  perturbed  ten 
mernbered  rings  with  the  perturbation  matrix  y  containing  only  one  non- 
-zero  off-diagonal  element,  ■=r\p^^  and  =  0  for  all  r  .  (For 

perturbation  to  Naphthalene  r  =  I  ,  8=6,  and  to  azulene  r  =  1  ,  s  =  5  , 

all  other  =  0  ).  Solution  of  the  equation  system  generated  by  equation 

(28)  gives  the  result 


Naphthalene  Azulene 


The  E(V)  vs  V  curves  iind  the  eigenvectors  can  be  obtained  as  before  but 
we  shall  give  only  the  curves  for  Naphthalene  which  indicate  how  complicated 
the  behaviour  can  be.  (Figure  6). 

Other  applications  of  this  type  have  been  made  of  the  method  but 
it  should  be  clear  from  these  examples  that  the  method  is  quite  useful.  We 
turn  next  to  illustrations  of  the  iterative  method. 


Figure  6.  Energy  levels  of  naphthalene  as  a  function  of  the  strength  of 
the  1-6  bond  perturbation. 
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IV.  APPLICATIONS  OF  THE  ITERATIVE  METHOD 

When  equation  (28)  of  section  II  yields  more  than  two  equations,  or 
perhaps  three  if  there  is  a  high  degree  of  symmetry  in  the  molecule,  it 
becomes  impractical  to  use  the  non-iterative  method.  We  shall  now 
illustrate  the  second  or  iterative  method  by  showing  how  the  molecules 
discussed  in  the  tables  compiled  by  Del  Re,  M^rtensson  and  Nordling 
can  be  treated. 


G.  Del  Re,  O.  Martensson  and  J.  Nordling,  Technical  Note  No.  29 
from  the  Uppsala  Quantum  Chemistry  Group  (1959). 


In  the  Hdckel  theory,  as  we  have  seen,  using  units  of  p  ,  we  can 

write 

(21)  G-H  -C  ’ 

and  we  know  how  to' obtain  this  inverse  by  the  method  of  Chebyshev  poly¬ 
nomials.  This  enables  us  to  solve,  by  iteration,  equation  (31)  of  section  II, 
1* 


(31)  (N)<c  (v-  V(&vX  =  0 

Because  of  the  cyclic  character  of  (G-  only  a  few  elements  of  the  inverse 
need  be  computed.  We  shall  list  below  the  quantities  needed  to  carry  out  the 
second  order  iteration  procedure  for  those  cases  considered  in  reference 
50.  These  are 


1)  (H)^  ^  on  the  basis  of  Hdckel  theory, 

2)  or  its  first  row  in  terms  of  Chebyshev  poly¬ 
nomials, 

3)  equation  (31). 

Although  these  equations,  when  written  out,  seem  rather  cumbersome 
the  actual  manipulation  of  them,  though  tedious,  is  not  difficult  and  the 
convergence  of  the  iterative  process  is  usually  quite  rapid  as  we  shall 
indicate  by  an  example.  In  the  figures,  the  large  dots  represent  large 
diagonal  perturbations  and  the  small  dots  represent  smaller  diagonal  per- 
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turbations  corresponding  to  the  inductive  effect  discussed  before.  In  addition 
to  these,  there  is  a  bond  perturbation  connecting  the  two  in  all  cases. 

It  should  be  obvious  from  the  information  shown  on  the  following 
pages  that  the  mathematical  structure  of  all  of  these  problems  and  many 
more  is  quite  similar.  This  is  especially  clear  from  the  form  of  the  A's 
and  B's  which  we  have  defined  as  an  aid  in  writing  the  equations  in  an 
orderly  form. 


(H), 


Cyclobutadiene 

3 


Perturbed  Butadiene 


0  10  1 
10  10 
0  10  1 
10  10 


(77) 


The  -1  in  the  and  positions  breaks  the  ring. 

*^11''  *^2^2  '  *  ^^13  "  ■  *^0^2  ^14 


=  'CiF; 


-1 


(78) 


AiCi  +  A2C2  +  =  0 

A2Cj+A4C2  +  A5C4  =  0  (79) 

A3C1  +  A5C2  +  A^C4  =  O 


^1 

»''u  -(Vu+V,2+  1)G„ 

-2(Vn^^l2-Vll)^12+2V^2G^3 

(80a) 

^2 

■  (^ll''^22  ^  12  ■  ^12^^12  ^22^13 

(80b) 

=  -1  +  V,,G,,  +  (V,2  -  1)G 

12 

(80c) 

^4 

,.V,2-(V%+vyG„- 

2^12^22^,2 

(80d) 

^5 

=  ^12^11  ^22^12 

(80e) 

^6 

"  -^11 

(80f) 
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Cyclopenta-dienyl 

radical 


I 


0  10  0  1 
10  10  0 
0  10  10 
0  0  10  1 
10  0  10 


'  V 


^11 

V12 

0  0 

-1 

^12 

V22 

0  0 

0 

0 

0 

0  0 

0 

0 

0 

0  0 

0 

-1 

0 

0  0 

0 

(81) 


—  (ri3  >  (^•y'*’ 


AiCi  +  A2C2  +  AjCg  =  0 
•^2^1  +  A^C2  +  A5C5  =  0 
AjCj  +  A5C2  +  A^Cg  =  0 


(82) 


(83) 


Note  that  although  the  expressions  for  the  A's  are  the  same  as  in  the  pre¬ 
ceding  case  the  matrix  elements  of  (S'  are  different  and  are  given  here 
in  terms  of  the  Chebyshev  polynomials  of  the  second  kind* 


Cyclopenta  -dienyl 


Perturbed  form  of  the  radical 


(H)e  same  as  in  the 
previous  case* 


V 


12 

0 

0 

12 


0 

0 

0 

0 

0 


0 

0 


0  0 
0  0 


1 


(84) 
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Note  that  Vjg  =  Vj2  ^55  “  ^ZZ  symmetry.  ^S"  has  been  given 

in  the  preceding  example 


B^Cj  +  B2C2  +  B3C5=0 

B2CJ  +  +  B5C5  =  0  (85) 

B3Cj  +  B5C2  +  B^C5=  0 


B,  =  V„-(V„  +  2V,2)G„-4V„V,2G,2-2v‘2G,3  (86a) 

=2“  ''l2-(''i,V,2+V,2V32)G„.(V„V22t2v22)G,2-V,2V22G,3  (86b) 

(86c) 


B3=  B2 


B. 


(V.o  +  V„)G,,  -  2V,,V,,G, 


-v"2°11-2Vi2V22Gi2-V22Gi3 


^6=  ®4 


(86d) 

(86e) 

(86f) 


Benzene 


pyridine 


i  -f 


0 

1 

0 

0 

0 

1 

^11 

^12 

0 

0 

0 

V12 

1 

0 

1 

0 

0 

0 

^12 

^22 

0 

0 

0 

0 

0 

1 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

> 

1 
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V22_ 

Gil - 

- 1,  f;'  ; 

&,3= -r.Fa-' 

Gn  = 

G-is-=  6-g  G-It=  6-ia 

(87) 

Bj^Cj  +  B2C2  +  830^  =  0 
®2^1  ^  ^4*^2  *  ° 

BjCj  +  B5C2  +  =  0 


(88) 
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G-i^r=-Tr'  j  G-ik  =  <ht.'= '  <  (,o) 


etc.  are  the  same  as  for  benzene 

P|C/i- P;^C{  rs.  0 

(91) 

C|  —  0 

> 

ill 

P 

(92a) 

d 

N) 

III 

(92b) 

°3  "  ^11 

■  ^ ii^n  ■  ^ii'^rr 

(92c) 

To  demonstrate  the  use  of  the  above  equations  and  to  indicate  the 
speed  of  convergence  of  the  method,  we  give  now  the  details  of  a  calculation 
on  a  perturbed  buta -diene -like  system. 

We  take  =  2  ,  V22  =  0.2  ,  V^2  “  *^“<1  obtain 

Aj  =  2  -  5.3Gjj  +  2.8  Gj2  +  0.6  Gj3 

A2=  0.3-0.66Gjj-0.19Gj2+0.2Gj3 

A3=  -1  +  2G^j  -0.7Gj2 

A4=  0.2-0.13G^^-0.12G,2 

Ag  =  0.3  Gjj,  +  0.2  Gj2 

■^6  ~  "  ^11 
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A  starting  value,  ,  of  E  is  chosen,  the  G's  obtained  and  the 
A's  calculated.  Any  two  of  the  equations  in  (79)  are  chosen,  Cj  is  set 
equal  to  unity  and  and  are  found.  We  obtained  the  following  re¬ 
sults  for  one  of  the  energy  levels  of  this  problem. 


TABLE  II.  Convergence  of  the  iteration  process  for  a  perturbed 
butadiene -like  system. 


Iteration 

E 

o 

Co 

•Is 

E*  (eq.  49) 

First 

3.0000 

1.0000 

1.0000 

2.77 

0.5359 

0.5361 

- 

0. 2009 

-0.2361 

0.0669 

Second 

2.77 

1.0000 

1.0000 

2.7798 

0.6031 

0.6033 

- 

0.2504 

0. 1046 

0.0905 

Third 

2.7798 

1.0000 

1.0000 

2.77996 

0.600021 

0.600021 

- 

0.247936 

0.089418 

0.089192 

Fourth 

2.77996 

1.00000 

0.999999 

2.779963 

0.599931' 

0.599958 

- 

0.247890 

0.089180 

0.089165 

We  see  from  this  table  that  the  convergence  is  quite  rapid  and  indeed  the 
second  iteration  would  suffice  for  most  purposes.  For  all  of  the  examples 
which  we  have  worked  out  the  second  or  third  iteration  has  given  us  the 
accuracy  we  needed.  Of  course,  the  number  of  iterations  required 
depends  on  the  accuracy  of  the  initial  energy  guess  E^  which  of  necessity 
should  be  chosen  different  from  the  unperturbed  eigenvalues,  but  in  most 
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case  8  one  has  a  rather  good  idea  of  what  this  should  be>  It  should  be 
mentioned  here  that  if  :my  of  the  eigenvectors  of  the  perturbed  problem 
has  C  ^  =  0  ,  then  that  eigenvalue  and  eigenvector  cannot  be  obtained 
by  this  second  order  iteration  procedure*  Since  one  is  dealing  with  finite 
matrices,  one  could  find  out  which  eigenvalues  are  lost  (excluding  those 
which  are  unperturbed)  and  apply  alternative  partitioning  procedure  to  get 
these  by  similar  procedures. 


V.  ,  DISCUSSION 

We  have  attempted  in  this  report,  to  indicate  how  the  methods 
described  in  reference  1  can  be  applied  to  problems  of  chemical  and 
physical  interest.  Here  we  wish  to  give  a  brief  discussion  of  the  relative 
merits  of  the  two  methods  and  to  point  out  some  possible  areas  for  applica¬ 
tion  of  them. 

First,  we  recognize  that  the  two  methods  are  basically  the  same; 
the  iterative  method  is  introduced  only  when  the  complexity  of  the 
perturbation  matrix  becomes  too  great  to  make  the  first  method 

tractable.  Hence  the  iterative  mbthod,  though  rather  complex,  is  much 
more  flexible  and  powerful  than  the  non-iterative  method.  The  usefulness 
of  the  non-iterative  method  can  be  extended  somewhat  by  going  to  first 
order  perturbation  theory  as  we  shall  now  indicate. 

Let  us  assume  that  the  Hamiltonian  matrix  of  our  problem  can  be 
written  as 

lHI=M,+v.+  v' 

where  QrQe  is  the  matrix  of  the  basic,  unperturbed  problem,  e.g.  benzene, 
is  a  large  perturbation  corresponding  to  one  of  the  problems  solved 
in  Section  III  and  is  a  small  perturbation  matrix  representing  diagonal 
and/or  off-diagonal  perturbations  on  the  substituted  ring  problem  whose 
matrix  is  M.tV.  P  then  from  first  order  perturbation  theory  we  have  that 

(94) 

We  have  already  solved  the  eigenvalue  problem  corresponding  to  Hat  VI 
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for  a  number  of  useful  cases  so  that  by  a  single  first  order  perturbation 
correction  we  can  extend  our  results  considerably. 

As  an  example,  we  consider  pyridine  with  both  inductive  effects 
and  bond  perturbations.  We  take  VJj—  =  0.  3  and  1 

in  the  perturbation  matrix  \D^  and  \^|  =1  in  the  matrix  with  all 

other  elements  in  both  matrices  equal  to  zero.  We  used  the  iterative 
method  to  find  out  the  exact  eigenvalue  to  the  accuracy  we  needed  and  in 
Table  III  we  compare  the  results  obtained  by  the  iterative  and  the  non- 
-iterative  perturbation  method.  Since  the  bond  perturbation  is  rather  large^ 
we  expect  that  most  of  the  error  in  ,  the  perturbation  value,  is 

introduced  by  it  and,  therefore,  we  also  show  the  results  of  the  two  methods 
when  0  but  s  0.  1  i.e.  inclusion  of  a  purely  inductive  effect. 


TABLiE  III.  Comparison  of  the  first  order  perturbation  results  with 
the  exact  results  obtained  by  the  iterative  methbd,  E^^for 
pyridine. 


Energy 

level 

j  O'l 

E^ 

El 

^1 

C 

2.  633 

2.6601 

2. 3125 

2. 3132 

B 

1. 370 

1.3550 

1.  3187 

1.3188 

A 

1.050 

1.0512 

1.0500 

1.0512 

A' 

-0.921 

-0.9057 

*0.6688 

-0.6684 

B' 

-0.950 

-0.9512 

-0.9500 

-0.9512 

C 

-1.982 

-2.0094 

-1.8625 

-1.8636 

It  is  obvious  from  this  table  that  an  inductive  effect  of  the  magnitude  con¬ 
sidered  here  can  be  included  by  first  order  perturbation  theory  with 
sufficient  accuracy  and  even  inclusion  of  the  bond  perturbation  in  this  way 
may  be  adequate  for  many  purposes.  In  any  case,  it  is  possible  to  con¬ 
sider  the  first  order  perturbation  energy  as  a  zeroth  approximation  in  the 
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iterative  method. 

The  first  order  perturbation  approach  should  be  quite  useful  in  all 
of  those  problems  where  there  are  a  number  of  large  diagonal  perturba¬ 
tions  which  are  almost  but  not  necessarily  exactly  equal  and  other  small 
diagonal  or  off-diagonal  perturbatiohs.  Many  compounds  and  series!  of 
compounds  of  biochemical  interest  should  fall  into  this  category. 

iThe  non-iterative  tnethod,  because  it  giyeS  a  closed  analytical  form 
for  the  E  =  E(V)  or  more  exactly  V  =  V(E)  .  relationship,  is  useful  for 
discussing  and  obtaining  insight  into  those  properties  of  simple,  cyolic 
systems  which  are  principally  determined  by  a  few  large,  localized 
perturbations.  The  details  introduced  by  small  perturbations  cam  then  be 
accounted  for  as  minor  corrections.  The  iterative  method,  on  the  other 
hand,  can  be  programmed  for  a  computer  and  its  value  will  be  particularly 
evident  when  a  series  of  calculations  on  related  systems  must  be  carried  out 
with  high  accuracy.  It  should  be  kept  in  mind  that  the  equation 
reduces  the  magnitude  d£  any  problem  considerably  because  the.  number  of 
non -vanishing  elemehts  in  \(?  becomes  the  measure  of  the  "size"  of  the 
problem. 

In  all  of  our  calculations  and  discussions  we  have  neglected  the 
effects  introduced  by  overlap  of  basic  atomic  orbital  set.  This  is  physically 
unrealistic.  The  second  order  iteration  procedure  can  be  extended  easily  for 
the  case  including  overlap.  Further  study  is  in  progress  on  this  aspect  of 
the  problem.  (See  Appendix  U). 

Turning  now  to  the  applications  of  these  methods  we  can  mention 
several  possibilities. 

a)  Further  applications  to  conjugated,  aromatic  systems. 

Besides  the  many  possibilities  which  exist  for  perturbation  on  the 
basic  benzene  ring,  there  are  similar  applications  to  five  and  seven  membered 
rings.  In  Section  lU,  we  showed  how  Naphthalene  and  azulene  can  he  treated 
as  perturbed  (crosslinked)  ten  membered  ring  and  we  have  treated  the  same 
molecules  with  substituents  at  selected  points.  It  should  eventually  be  possible 
to  treat  many  more  cross-linked  chains  by  these  methods. 

The  application  of  the  method  to  molecules  of  biochemical  inport.ance 
should  be  quite  profitable.  Even  the  very  simple  quantities  obtained  for 
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these  molecules  by  Hflckel  theory  such  as  the  energies  of  the  highest 
occupied  molecular  orbital  (HOMO)  and  the  lowest  empty  molecular 
orbital  (LEMO)  have  been  of  great  value  in  obtaining  insight  into  the  elec¬ 
tronic  nature  of  biochemical  processes.  Here  one  is  not  so  much  inter¬ 
ested  in  the  exact  energy  values  but  rather  in  the  variation  of  the  energies 

51) 

among  a  series  of  molecules  of  the  same  general  structure  It  is 
also  important  to  realize  that  with  these  methods  it  is  not  necessary  to 
calculate  all  of  the  energy  levels  simultaneously  as  in  the  case  with  the 
usual  solution  of  the  secular  equation  of  the  problem.  One  can  concentrate 
attention  on  the  energies  of  the  HOMO  and  the  LEMO  and  follow  their  varia¬ 
tion  with  molecular  structure  and  forget  about  the  other  energy  levels. 


51) 


B.  Pullman  and  A.  Pullman,  Revs.  Modern  Phys.  32,  428  (I960) 


b)  One  dimensional  problems. 

The  use  of  one  dimensional  models  as  prototypes  for  three  dimen¬ 
sional  problems  is  very  common  and  valuable  in  physics  and  chemistry 
even  though  three  dimensional  problems  may  exhibit  effects  not  occurring 
in  the  one  dimensional  analogues.  We  consider  it  quite  possible  that  such 
problems  as  the  electronic  and  vibrational  properties  of  a  linear 
chain  with  impurities,  one -dimensional  alloys,  anti -ferromagnetic  linear 
chain  and  other  one -dimensional  order -disorder  phenomenon  can  be  treated 
either  exactly  or  in  good  approximation  with  these  methods. 


52) 


See  reference  42. 


53) 


See  reference  43. 


c)  Three  dimensional  problems. 

The  extension  of  the  one -dimensional  calculations  to  two  and  three 
dimensions  is  not  trivial  although  one  has  an  intuitive  feeling  that  the  cyclic 
symmetry  which  permits  the  inverse  of  the  corresponding  matrix  to  be 
found  by  the  method  of  Chebyshev  polynomials  can  be  displayed  in  three 
dimensional  problems.  This  will  probably  be  easy  in  some  cases,  difficult 
in  others.  In  any  case,  we  can  be  certain  that  in  solid  state  physics,  even 
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the  simple  calculations  which  are  based  on  assumptions  similar  to  those 
made  in  the  Hdckel  theory  of  molecules  will  be  needed  for  some  time  to 
come  and  hence  we  need  not  be  embarrassed  by  the  crudity  of  our  assump¬ 
tions. 

We  have  work  underway  on  some  of  these  applications  and  we  shall 
try  to  report  further  on  them  at  a  later  date. 
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APPENDIX  I 


CHEBYSHEV  POLYNOMIALS 


The  Chebyshev  polynomials  of  the  first  kind  are  defined  by 
=  Z  Cos  nb  ^  Z  CoS  (.rt  om. 

where  S,Co5  6. 

The  polynomials  satisfy  the  following  recurrence  relation 


The  polynomials  of  any  arbitrary  degree  n  can  be  generated  by  the  above 
recurrence  relation  with  the  starting  values  . 

The  Chebyshev  polynomials  of  the  second  kind  are  defined  by 

=.  Sen  O^-t-QQ  ^  SCm.  [7?v+-fj 
SC-n& 


where  x  =  Z<u>sb . 

The  polynomials  'C-riOQ  satisfy  the  recurrence  relation 
-  x^»v-(  Cx;)  -  ^n-2.Cx-7 


with 


and 


Cx)  =1  X  . 


The  relation  between  is  given  by 

'^n  (*•)  ^  —  %-n-^ 

The  explicit  representations  of  and  in  powers  of 

'V.W=  X--(V)  . 

=  X-  -cxrr') 


X  are 
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APPENDIX  ri 

INCLUSION  OF  OVERLAP 


We  have 

^  E- ^ 'll  So’+ sj) ^  (Ai) 

where  §  is  the  perturbation  matrix  on  which  is  the  overlap  matrix 

characteristic  of  the  cyclic  unperturbed  problem, 

(A2) 


We  can  write  equation  (Al)  as 

C  =.  I^E('1l-hS0  -  OQ  (V-ES)<E 


(A3) 


Drita.  Wfej  =  V’-ES  ;  K  - 

and  note  that  tvP  =: 


SO 


(C=  IKW>C 


(■D  -  ll<;w?)c  -  0 


(Ad) 


The  matrix 


IK 


can  be  written  as 


K=  He.  +'S,nii+stni+ 


t  -  IHP 


where  Sj  is  the  overlap  integral  between  nearest  neighbours,  between  second 
nearest  neighbours  and  so  on.  So  the  matrix  can  be  inverted  in  terms 

of  Chebyshev  polynomials  (Section  III  of  reference  44)  and  we  can  apply  the 
non-iterative  method  or  iterative  method  as  before  with  slight  modification,  • 


To  apply  partitioning  technique,  we  convert  the  matrix 
into  a  hermetian  matrix  by  premultiplying  by  .  So  we  get 


ffl-IKW) 


(W  -W  kM)€  =  0 


(A5) 
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Defin*  (M)  =  (W  -  UP  fK  ^ 

(MIC  =0 


(A6) 


Following  the  procedure  outlined  earlier,  we  construct  a  second  order  itera¬ 
tion  scheme 


E„ 


0^, 


where  the  prime  denotes  differentiation  with  respect  to  £  .  Now 

M=  (V-ES)  -(V’-ES)  K  (V-ES) 

-5 +SK(V’-ES) -t-CV’-ES^IICS  -  (V-ES>k'CV»-ES) 


Now 


CeC'i+so-wrj’ 

!K'=  -1X('1/+s.)IK 


So 


ITD^'  -5  +  S(K(V-ES)-Ky-ES3KS  +  {V-ES)4!C(/PtS^IXCV':E5) 

For  exact  E  euid  O  ,  we  have  c=  •  But  for  arbitrary 

E„,  (P-e,  we  have  in  general 

do  =  (Ko^vPoC  ; 


-e;^s<Co +cc^sci.>> +d^sc.  + 

Defining  Sff. = £  ^  finally  we  get 
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So  the  final  iteration  scheme  becomes 

^  ^  ^  ,  _ _ _ _ _ 


Following  reference  1,  we  can  also  consider  the  non-iterative  scheme. 
For  one  diagonal  element  perturbed,  we  get  from  the  eqviation 

cU  [4-lKW]  =  o  ^ 

I  —  ktt  wii 

which  yields  the  energy  levels  E  in  terms  of  V^j,  .  In  a  similar  way,  we 
can  derive  an  analogous  formula  for  one  bond  being  perturbed  only.  These 
equations  are  slightly  more  complicated  than  those  in  the  case  of  overlap 
neglected  but  still  simple  enough  to  be  quite  useful. 
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